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THE ITERATION OF FUNCTIONS OF ONE VARIABLE. 

By Albert A. Bennett. 

A large part of the theory of the iteration of functions of several variables 
is illustrated by the relatively simple case of functions of one variable, and 
it is to the study of this case that the present paper is devoted. Several 
distinct classes of topics immediately suggest themselves in this connection. 
One of these classes comprises subjects of a purely formal nature, and is 
concerned only with expansions in the neighborhood of certain significant 
points. Another class of topics deals with the iteration of a real function 
within an interval. Still a third treats of the total behavior of the trans- 
forming function B{x), when the given function A\(x) is analytic and 
defined for all values of x. Except for certain simple cases the function 
B(x) cannot be one-valued. The present paper deals almost exclusively 
with the first two of the three classes of subjects just mentioned. In the 
formal problem, certain new points of view, and new results and formulae 
are obtained, although a large proportion of this part is a systematization 
of results already known but largely unrelated. The second part of this 
paper deals with the iteration of a real function, a question which seems to 
have been neglected up to date.* The third topic mentioned above, has 
not here been touched upon, except in the case of certain simple rational 
functions. A few known theorems concerning the iteration of functions 
of a single variable, have been omitted, in view of the fact that they seem 
to have little significance in connection with the subject of the iteration of 
functions of several variables. Only those subjects have here been dis- 
cussed which appear as natural introductory material to the more general 
case. 

For references to those parts of the formal problem which have already 
been treated in the literature, one may consult the following : 

Pincherle, on " Functional Equations and Operations," in the Ency- 
klopadie d. Math. Wiss., II, All, and Encyclopedie d. Sci. Math., II, 26. 

E. Schroeder, Math. Ann., 2 (1870), p. 317, and 3 (1871), p. 296. 

J. Farkas, Journ. de Math. (3), 10 (1884), p. 102. 

G. Koenigs, Bull. Sc. Math. (2), 7 (1883), p. 340, and (3) 1 (1884), 
suppl., p. 14. 

* Since this paper was written, there has been a report on one phase of the case of real 
iteration by Mr. J. F. Ritt. Cf. Am. Math. Soc. Bull., 21 (1915), p. 379. 
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24 ALBERT A. BENNETT. 

E. Podetti, Giorn. Mat. (2), 4 (1897), p. 264. 

C. Formenti, Reale 1st. Lomb. Rend. (2), 8 (1875), p. 275. 

A. N. Korkine, Bull. Sc. Math. (2), 6 (1882), p. 228. 

C. Bourlet, Toulouse Ann. (1), 12 (1898), mem. no. 3, p. 1-3. 

A. Grevy, Thesis, Paris, 1894, and Ann. Ec. Norm. (3), 11 (1894), p. 287. 

L. Leau, Thesis, Paris, 1897. 

E. M. Lemeray, C. R., 125 (1897), p. 524, also S. M. F. Bull., 26 (1898), 
p. 10. 

0. Spiess, Math. Ann., 62 (1906), p. 226. 

M. Koppe, Die Iteration des Sinus . . ., Berlin, 1909. 

Pincherle and Amaldi, Le Operazione Distributive, Chap. XIV. This 
may be consulted in connection with part of the use here made of matrices. 

Note. The matrices and iteration discussed in this paper have no rela- 
tion to the matrices and " Schapira's Iteration " treated in several recent 
papers by L. v. David, in the J. fur Math. 

PART ONE. THE ITERATION OF POWER SERIES. 
Definition and Examples of Iteration. 

2. Let there be given an analytic function, A\{x). The function 
Ai[Ai(x)], we shall define as A 2 (x), and in general for n, a positive integer, 
we shall define A„(x) as identical with An-i[Ai(x)], or, what amounts to the 
same thing, as identical with Ai[A n -i(x)]. We may similarly denote the 
inverse of Ai(x) by A-\{x). For two given functions Ai(x), and Fi(x) 
we may define a (?i(x), so that Gi(x) = F_i{Ai[Fi(x)]}, or, in other words, 
so that Gi(x) is the transform of Ai(x) through F-i(x). The Gi(x) and 
Ai(x) will then be so related that for any positive integer, or even negative 
integer value of n, G n {x) = F-i{A n [Fi(x)]}. In particular, we may 
replace the consideration of y = A\(x) by that of rj = (?i(£), where 77 is 
obtained from y, and £ from x, by the same linear fractional transformation. 

It is natural to inquire whether we might not suppose Ai(x) and Gi(x), 
as given, and determine a transforming function F, such that Gi(x) 3= 
F-i{Ai[Fi(x)]\, or, as it may be written, Ai[Fi(x)] = Fi[Gi(x)]. In 
particular, we may examine the case in which Gi(x) is simply x + 1. The 
function F associated with A by this particular choice of G, we shall call B. 
We seek therefore to determine a function B, such that 

(1) AJLB(x)] = B(x + 1). 

We might equally well have started with Gi(x) as ex, where c is a constant, 
=t= 0, 4= 1, and have sought a function E(x), such that 

(2) A 1 [E(x)] m E(cx). 
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The two problems are not, of course, independent. We need merely replace 
x in (2) by c x , to obtain 

A 1 [E(c*)] B E(c*+ l ), 

so that E{c x ) is of the form of B(x). If a function B{x) be found, we may 
define A n (x), as identical with Bi[n + B_i(x)], for n not only an integer, 
but, indeed, for all complex values of n, and similarly for E(x). The 
definition of A n (x), will depend to a considerable degree upon the particular 
function B or E, that is selected, except for n an integer. We shall define 
A n {x) as the nth iterate of Ai(x), and the process of repeatedly substituting 
Ai(x) in place of x, or, in general, of determining A n (x) from Ai(x), as the 
process of iteration. 

We shall now give a few simple examples of iteration, for which no 
extended discussion will be necessary. For the degenerate case of Ai(x) 
equal to a constant, equation (1) defines B(x) as also equal to a constant. 
For Ai(x) = x, B{x) may be any periodic function P(x), with unity for a 
period. For Ai(x) = x + 1, B(x) is of the form x + P(x), where P{x) 
is defined as above. For Ai(x) = cx, + c # 1, .B(x) is of the form 
gCx+pomog^ jd^-j ag b e f ore> This includes such expressions as ac x , where 
o is an arbitrary nonzero constant, and more generally, it includes an exten- 
sive class of elliptic theta functions. For Ai(x), an arbitrary linear frac- 
tional function, we obtain, either directly, or by reducing to a normal 
form, a class of associated functions, B(x), containing an arbirtary periodic 
function P(x), with unity for a period. 

If Ai(x) be not a linear fractional function, then either Ai(x) itself, or 
else its inverse, A_i(x), must be multiple- valued, which would suggest that 
no one- valued function, B(x), could exist for such a function A x (x). This is 
not, however, the case. For example, if we define Ai(x) by the implicit 
equation: 

Coo + CioX + C 20 X 2 



+ Coi^.i(x) + c n xAi(x) + Cnx-Axix) 
+ c 02 Ai 2 (x) + c u xAi 2 (x) + c 2 2X 2 A 1 2 (x) 



= 0, 



where d, = Cji, i, j, = 1, 2, 3, and where the determinant of the coefficients 
is different from zero, we shall have a non-degenerate quadri-quadric rela- 
tion between Ai(x) and x. For this choice of Ai(x), there exists among the 
solutions of (1), a B(x), which is a single-valued elliptic function. The 
corresponding A n (x) satisfies an analogous quadri-quadric relation for 
every particular choice of n, whether real or complex. This is the problem 
which arises in connection with Poncelet Polygons, and the addition 
theorem of elliptic functions.* 

* Cf. Cayley, Elem. Treat, on Ell. Funct., p. 340, and Halphen, Traite des Fonc. Ell., 
vol. II, Chaps. 9 and 10, and an article by the author in Annals of Math. (2), 16 (1915), p. 97-118. 
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The question arises as to what is the most general algebraic expression 
Ai(x), which is defined by equating to zero a polynomial in A n (x) and x, 
whose coefficients vary with n, but whose degree is, in general, the same, 
and for which there exists a B(x), satisfying (1), and one-valued all over the 
complex plane. This problem involves the question as to what types of 
algebraic curves admit of one to one continuous transformations into them- 
selves. As is well known, the rational and the elliptic curves are the only 
ones admitting such transformations. The problem just proposed reduces 
to the study of addition theorems of rational and of elliptic functions. 
We shall not go further into this matter, but shall proceed immediately 
to a more general problem. 

Significance of the Formal Iteration of Power Series. 

3. If we confine ourselves to the immediate neighborhood of a finite 
value of x, the equation (1) imposes no restrictions on B(x), and is without 
significance. It is only when x is in the neighborhood of infinity, that we 
can expect to draw formal inferences with regard to B(x), as a consequence 
of the behavior of Ai(x). We shall seek to determine a solution B{x), 
which at least admits asymptotic representation along the positive real 
axis, using the term "asymptotic representation," in its widest sense. 
In particular, we shall require that lim B(x) exists. By replacing Ai(x), 

if necessary, by a linear transform of it, we may secure that lim B(x) = 0, 

or from (1), that lim Ai(x) = 0, when x approaches zero along a suitable 

path. There is another reason, also, which induces us to require that 
lim Ai(x) = for at least one method of approach, a reason, which is 

not explicitly related to the nature of the transforming function B(x). 
The problem of computing the coefficients of a power series representation 
of Az(x) from the coefficients of a series representing Ai(x) requires that 
at least one point (e) be known, either finite or infinite, for which A(e) = e, 
unless we are to have recourse to infinite summations which cannot be 
convergent except when the circle of convergence for the series Ai(x), 
includes the point corresponding to the constant term of this series. 

The derivation of B(x) from Ai(x) involves two steps; the first, that of 
the purely formal derivation of the coefficients, the second, that of proving 
the convergence or divergence of the series formally obtained for B(x). 
The first of these steps is not, however, entirely devoid of significance, even 
when the resulting series is divergent. The series 

A^x) 33 a n x + a 12 x 2 + ■ ■ ■ + a ln x n + ■ •• 

may be divergent, and yet it may well happen that there exists an analytic 
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function F{x), with an isolated essential signularity at the origin, for which 

along a suitably chosen path of approach. This is the case presented by 
the classical illustration of y = e~ llx \ where the path may be taken as either 
real semi-axis. Were we to choose infinity instead of the origin as the 
point approached, we should have the usual case of asymptotic convergence.* 
The set of coefficients an, au, ai 3 , • • •, ai„, • • •, determines what may be 
called a differential element of infinite order. Two distinct functions may 
admit of the same asymptotic representation for a given method of ap- 
proaching infinity, or to use more general ideas, they may have the same 
differential element at a given point, for a given path of approach. The 
two functions then represent curves which along a given path in the complex 
plane, have contact at the given point, of infinite order. So far as the 
formal determination of the series is concerned, we may suppose ourselves, 
throughout, as interested only in "asymptotic" forms, i. e., in differential 
elements for which questions of convergence are without significance.* 
It has been shown and it may be immediately verified that the analytic 
functions represented by 

where 1 < b n > | a n |, b n and a„ being real, have in general, an essential 
singularity at the origin, but have for real values the common differential 
element (a h a 2 , • ■ •, a n , • ■ •), at the origin. 

It would, of course, be possible to consider series Ai(x), with no absolute 
term, but with terms containing other than merely integral powers of x. 
These, however, we shall not investigate. Even in the case of direct power 
series, a classification of types will be found necessary. We might enunciate 
an apparently arbitrary classification, justified a posteriori by the different 
conclusions obtained. We shall however find that the types which arise 
are essentially those which are suggested by regarding iteration as a case of 
collineation in an infinite number of dimensions. It is this point of view 
that we shall now present. 

The Iteration of Power Series as Matrix Multiplication. Classification 

of Types. 

4. For many purposes, it is more convenient to replace the nonhomo- 
geneous relation y = Ai(x) by the pair of homogeneous equations, obtained 

* Cf. W. B. Ford, S. M. F. Bull, 39 (1911), p. 347. 

tCf. Kasner, "Conformal Geometry," Fifth Int. Cong, of Math. Proc, vol. II, p. 81 ff. 
The term, "differential element of infinite order," has been discussed by Kasner. See also 
J. F. Ritt, Bull. Am. Math. Soc, 21 (1915), p. 379. 
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by replacing x by Xi/x , and y by yi/y - We shall indeed write 

yi = XoAl (S)' 

Uo = x . 

We define Ai (m) (x) as the polynomial of the rath degree in x, which is ob- 
tained by neglecting all but the first to terms of Ai{x). Now the first to 
terms of A n (x), for n a positive or negative integer, depend only upon the 
first to terms of Ai(x), that is, A n {m) {x) is determined from Ai(x) by using 
only Ai (m) (x). 

Every binary form of the wth degree in y , y\, will be itself a binary form 
of the wth degree in xo, Xi, together with an infinite series, every term of 
which contains x to some positive power in the denominator. We shall 
now define £i (m) and ??, (m) as follows: 

£,-<"»> = av-xi" 1 - 1 ' i = 0, 1, 2, • • •, to, 
i/, (m) = integral part of series for i/o i -2/i m-t -= Xo i [x ^4.i(xi/a;o)] ,n_< . 

Here £, (m) is a product transform of the variables x , and Xi, and rii im) is 
related to the product transform of y and y\* The expressions ij< (m) , 
i = 0, 1, 2, • • •, to, are linear in the to + 1 homogeneous variables, £, (m) , 
i = 0, 1, 2, • • •, to, and in particular r? (m) = £ (n) . The matrix of this 
linear transformation from the (to + l)£ (m) 's to the (to + 1), i7 (m) 's we may 
represent by Ai> m) . 

For any other series y = H(x) we may make an analogous discus- 
sion, and denote by i/ (m) the analogous matrix. The first to terms of 
z = H[Ai(x)], where H(0) = 0, depend upon Ai (m) and H<- m) alone. Let 
us now define U im) as the integral part of the series for zo'.zi" 1-1 = 
yo^yoHiyi/yo)]" 1 -', where zi = y H(yi/y ), and z = 2/o- We may deter- 
mine the first to terms of z = H[Ai(x)] in three steps, by first taking the 
to leading terms of y = A\{x), secondly, taking the to leading terms of 
z = H(y), and thirdly, substituting the first to terms of Ai(x) for y in the 
first to terms of z = H(y), and taking the to leading terms of the result. 
This corresponds to considering, first, the symbolic equation ?? (m) = .Ai (m) £ (m) , 
secondly the symbolic equation f (m) = H < - m) ri (m) , and finally, substituting 
•Ai (m) £ (m) for 77 (m) in the equation for f (m) . Since Ai<- m) and // (m) are the 
matrices of linear transformations, the result of the substitution yields 
merely a new matrix of the same type, which is the usual matrix product, 
/f (m) Ai (m) . If in particular // (m > = A n -^ m \ then we see that A n (m) = 
An-i (m) Ai<- m) , where the product is the usual matrix product. The signi- 
ficance of this relation lies of course in the fact that the initial first minor 
of i4„ (m+1) is simply i4„ (m) itself. 

* Cf. A. Htirwitz, Math. Ann., 45 (1894), p. 388, for the question of product transforms. 
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The relation Gi(x) = F-i{Ai[Fi(x)}] is equivalent to the equation 
Gi (1B > = Fi (m) '- 1 i4i (m) Fi (m) , where F^-- 1 is the inverse of Fi (m) in the 
notation of matrices. This is the type of equation that arises in the theory 
of the classification of the linear transformations Ai (m) . Although the 
matrix Fi (m) is not of the most general form possible, still it is not surprising 
that many of the theorems concerning the normal forms of linear trans- 
formations have here an application. For the sake of concreteness, we 
now give the form of i4i (5) 

10 

an aw ai3 a u au 

an 2 2onOi2 2onai3 + au 2 2auau + 2auai 3 

an 3 3a n 2 ai2 3an 2 ai 3 + 3anai 2 2 

an 4 4au 3 ai 2 

an 5 

The case of Ai (5) is typical of A n (m) to the extent that the first row contains 
only zero except for a one in the initial position, the elements below the 
main diagonal are all zero, while the elements in the main diagonal are the 
successive powers of a„i. The roots of the characteristic equation of the 
matrix Ai (m) are clearly flu', i = 0, 1, 2, • • •, m. We shall now examine the 
normal forms for the linear transformations whose matrices are Ai^K 

The first case to consider is that in which all of the roots of the charac- 
teristic equation are distinct. This is the case in which a n is not an integral 
root of unity. We may without exception in this case* find matrices 
£(»») by the method of undetermined coefficients, for which 



Ai<« = 



A^ m ~> = E {m) <- l A^ m) E { 



m) 



where Ai<- m) has the same main diagonal as Ai im) , but has only zeros above 
as well as below this diagonal. The matrix £ (m) has a single arbitrary 
constant which is determined completely, if, for example, we put e\ = 1. 
The only invariant of i4i (m) as m takes on successively greater integer 
values, is in this case an. 

We next come to those nonsingular cases in which all of the roots of 
the characteristic equation coincide. This occurs when on = 1. We may 
now have as many as two algebraic and one arithmetic invariant. It is 
impossible to have more, however, no matter how great m be taken. When 

* This first case has been treated by nearly all of the writers on the subject, Schroeder, 
Korkine, Farkas, Koenigs, etc., already referred to. This use of matrices A„ im) was made 
for the first time by H. v. Koch, Bil. t. Sv. Vet. Ak. Hand. 1, Math. 25 (1900), m<*m. no. 5, 
p. 1-24. 
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an = 0, for every i > 1, we have the case of the identity. We shall have 
in all other cases, an arithmetic invariant k, such that au = 0, 1 < i ^ k, 
while ai, fc+i 4= 0. This is an invariant in the sense that no matter how F (m) 
be chosen, the positive integer k for d (m) where Gi (m) = F (m)>_1 i4i (m) F (m) , 
will be the same as for .Ai (m) . By the use of undetermined coefficients we 
may find a matrix £ (m) such that 

where this time we mean by Ai im) the matrix obtained from a series of the 
form 

Ai(x) = x + a uk+1 x k+1 + ai,2k+ix u+1 + ■ ■ ■ + 5i, ab+i***- 1 + ■ • •, 

and where (l/x)Ai(x) is a power series in x k . Both oi,»+i and oi, 2*4-1, are 
invariants of Ai(x) in the sense that any two series A\{x) equivalent to the 
same series Ai(x), must have the same oi,*+i and 5i ,2*4-1.* 

Of the nonsingular cases there remain to be considered those in which 
some but not all of the roots of the characteristic equation coincide. For 
these cases an must be an integral root of unity and we may write an* = 1, 
where h is the smallest positive integer for which this is true. But if Ai(x) 
be of this form then Ah(x) will be of the type already mentioned, since it 
will commence with x + • • • . Thus we have as invariants at least two 
integers h and k, and two algebraic invariants a,h,k+i and Si, 2*4-1. Here k 
must be a multiple of h but need not be h itself. Furthermore by the use 
of the method of undetermined coefficients, these are found to be the only 
invariants. 

The only other cases which can arise are the singular cases, and in these, 
an = 0. If every succeeding an, i = 2, 3, • • • , vanishes, the expression 
Ai(x) vanishes identically and is independent of x. Let us suppose then 
that Oii = 0, i = 1, 2, • ■ •, k — 1, while Oi* ={= 0. In this case k will be 
an arithmetic invariant, and this is the only invariant, f When once a 
normal matrix A^ m) has been selected, the remaining problem is to express, 
if possible, A n <- m) as a matrix whose elements are explicit analytic functions 
of n. 

Synopsis of Types. 

6. On the basis of the above discussion we are prepared to enunciatej 
a classification of power series of the form 

Ai(x) = Onx + ai 2 a; 2 + • • • + OirX r + • • • 

from the point of view of iteration, into eleven types. We shall later 
examine the various cases in greater detail. 

* For this and the following case, see Leau, loc. cit. Cf. also Kasner, loc. cit. 
t This is the case particularly investigated by Grevy, loc. cit. 
% This is given here for the first time. 
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Case I. an is not an integral root of unity and does not vanish. 

Type la, | a n | < 1, 4= 0. 

Type 16, | a n | > 1. 

Type Ic,* | an | = 1, but a n h 4= 1, h = 1, 2, • • •. 

Case II. a u = 1. 

Type Ha, au = 0, t = 2, 3, 

Type 116, an = 0, i = 2, 3, • • •, k, but ai, k +i 4= 0, and 2ai, 2k +i = 
■ (k + l)ai 2 ,it + i, where ai, k+ i and ai, 2 *+i are the invariants of the normal 
form already mentioned. 

Type lie, an = 0, i = 2, 3, •••, k, but ai,* + i 4= 0, and 2ai, 2 *;+i =1= 
(A; + l)5i 2 , jfe+i, where 5i, i+i and 5i, 2*+i are defined as above. 

Case 7/7. an 4= 1, but On' = 1, where /i is the smallest positive integer 
for which this is true. 

Type Ilia, a*; = 0, i = 2, 3, 

Type III6, a hi = 0, i = 2, 3, ••-,&, but a,,, fc+ i 4= 0, and 2a h , 2 k+i = 
(k + l)ah 2 , k+i, these being defined as above. 

Type IIIc, a hi — 0, i = 2, 3, • • •, k, but a h , k +i 4= 0, and 2a fc|2 ifc+i 4= 
(A; + l)a A 2 ,*+i, defined as above. 

Case IV. a n = 0. 

Type IVa, a u = 0, i = 1, 2, 

Type IV6, an = 0, i = 1, 2, • • •, k - 1, but a lk 4= 0. 

Commutative Series. 

6. If the approximate structure of a series for A n (x) be once found in a 
more or less indefinite fashion, it is sometimes possible to obtain the coeffi- 
cients of A n (x), as explicit functions of n, without first reducing A x (x) to 
a normal form Ai(x). This is done by using the following considerations. 
If a series A n (x) exists, it must satisfy the commutative condition 

AnlAiix)] = Ai[A n (x)]. 

Whenever the solutions C(x), of the equation, 

(3) C[Ai(x)] « MC(x)] 

form a one-parameter family, we may seek to identify this parameter with 
the n of A„(x). If, however, the C's form a discrete set, or contain more 
than one parameter, we cannot obtain A n (x) by this means. In terms of 
matrices, we have the equation, 

Q(.m)J^(m.) _ ^J (m)^(m) 

which obviously greatly restricts C (m) . It is found that except in those 
cases in which an is an integral root of unity, there always exists a unique 

* For questions of convergence this yields three quite distinct subcases, as we shall see later. 
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one-parameter family of series C(x), commutative with Ai(x), in the above 
sense, provided that in the cases for which an = 1, we require also that 
Ci = 1. In the cases in which an = 0, the series C, although existing, 
are not, in general, in integral powers of x. 

The identification of the parameter, with w, is, of course, not unique, 
but involves an arbitrary periodic function of period unity. A particu- 
larly simple choice is in each case almost immediately obvious. For 
instance, to illustrate the most common case, let 4= | «n | 4= 1, and let 

C(x) = C\X + CiX 1 + c s x 3 + • • • 

while A\(x) is denoted as heretofore. Then C[A\(x)] will be of the form, 

CidnX + C1CI12X 2 + CidizX 3 + • • • 

+ c 2 an 2 x- + 2c 2 anai2X i + ■ ■ • 

+ Csdi^X 3 + • • ■ 

+ ••• 

and Ai[C(x)] will be obtained from this by interchanging a u and c,-, i = 1, 
2, 3, • • • . Equating coefficients of like powers of x, we obtain, 

Cidi = OuCi, 

Cittu + C 2 On 2 = OnC 2 + «l2Ci 2 , 

CiOis + 2c 2 anai2 + c s a n 3 = onCs + 2a n CiC 3 + ai 3 Ci 3 , 



Hence 



Ci (ci - 1) 
an (an — 1) 



aT; (an 2 - 1) ° 13 + 2 «7! (aTT^l) (Cl " an)ai2 • 



If C(x) is to be identified with A n {x), we must have for integer values of n, 
ci = an". The most general form of Ci is then, an n+P(n) , where Pin) is a 
periodic function of period unity, vanishing for n an integer. The most 
obvious choice of Ci is undoubtedly that in which P{n) vanishes identically, 
so that for all values of n, c x = an". Every succeeding c iy i = 2, 3, • • • , 
is then completely determined. 

Discussion of the Different Types. 

7. For any type under Case I, we may always determine a series E(x), 
and if we put e x = 1, it will be determined uniquely and in such a way, 
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that 

A 1 (x) m E-^auEix)]. 

By writing E-i(a n x ) = B(x), we obtain a B(x), for which 

A n (x) B JB[» + B_i(x)], 

as was desired. This gives A n (x) as an explicit function of n. We may 
obtain, however, another form of A„(x), which proves for many purposes 
much more convenient than that just given. It is to be expected that 
by the use of the expressions Ai(x), A 2 (x), • • ■, A p (x), • • •, where p is a 
positive integer, we may obtain a sort of interpolation formula, giving 
A n (x) for all values of n. We shall seek indeed a series of the following form 

A„(x) = J2cnUi(x), 

i=0 

where C/»(x) is independent of n, and is a power series in x commencing 
with a term in x m , expressible as a linear function of A (x), Ai(x), • • •, 
Ai(x), with coefficients which are themselves polynomials in an with 
numerical coefficients, and where a» is a function of On and n alone. These 
conditions serve to determine the U's and a's as follows: 

U (x) = A (x) s= x, 

£/i(x) = Ai(x) - anA'o(x), 

f/ 2 (x) = [A 2 (x) - 011^1(0;)] - an^A^x) - anAo(x)], 



U,(x) = C/ i _i[A 1 (x)] - On^-iW. 

The coefficients, a,-, satisfy the recursion formula, 

1 (a u "- m - l)ai-i 
a «- = a7x (a u ' - 1) ' a ° = au ' 
Then we will have 

A.ix) . au- [ £7„(x) + i ('"'r^ p " W 

(4) +Afa^ ( ?"T",y g»W+-l- 

On 2 (a u - 1) (an 2 - 1) ' w J 

To prove that the A„(x) determined by (4) is actually identical with that 
determined by the use of the series B{x), we need only note the following 
facts: 

1. Since Ai(x) commences with the term o u % the series C/,(x) must 
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necessarily contain no power of x less than the (i + l)st, as is readily 
verified by the use of the recursion formula for the U's. 

2. The equation (4) is an identity for n, a positive integer, say i — 1. 
For then, 

1 (a u ° - !)«<_! 



Oii = 



an (an*' - 1) 



will vanish, and all succeeding a's, while the non-vanishing terms add up 
to exactly Ai(x), as is readily verified. 

3. As a consequence of the two facts just mentioned, the coefficient of 
x p , in A n (x), p a positive integer, and n arbitrary, will depend upon only 
the first p, U's. It will therefore be a polynomial in an" and contain n 
in no other way, the coefficients being rational in an and integral, rational 
in the first p coefficients of A x {x). 

4. The coefficients of x p in A„(x), as obtained by the use of B(x), will 
also satisfy the conditions just stated, and since these two polynomials in 
an" must coincide for an infinite number of values of an", viz. for n = p + 1, 
p + 2, • • • , the polynomials are identical. Thus the two determinations of 
A„(x) coincide. 

Another* and somewhat analogous representation for A n (x) may be 
obtained from (4), in which the variable n occurs in a more convenient form, 
but in which the coefficient of x p , p > 1, is in general an infinite series in an. 
This is defined in terms of V's, as follows: 

V (x) = A (x) = x, 

Vi(x) = Ai(x) - anAo(x), 

V 2 (x) = [A 2 (x) - anAiix)] - an[Ai(x) - a n A (x)] 

= A 2 (x) - 2a 1 i4i0) + a n 2 A (x), 



(5) Vi(x) s V^A.ix)] - anFi-i(x) 

= Ai(x) - (I JauA^ix) + f * Jan 2 ^.»-2(a;) - • • • * a^A^x), 

A n (x) = au n [V (x) + (^±V 1 (x) + (l^ 2 V 2 (.x) + ... 
In Type Ha, Ai{x) is the identity. This case requires no further men- 



Equations (4) and (5) are given here for the first time. 
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tion. Type life is a well-known form discussed by many of the authors 
already mentioned. By the use, for example, of the method of undeter- 
mined coefficients, the Ai(x) may here be transformed through a suitable 
F(x) into a series of the form x times a power series in x k , the leading term 
being exactly x. This applies also to Type He. There is a wide degree of 
freedom in the choice of the series. In fact only the coefficients of the terms 
in x, x k+1 , and x ik+1 are preassigned. For Type lib, we may in particular 
take the series to be in the form 

y = Ai(x) = 



a/1 — kai, k +ix k 
If in this expression we replace y by y llk , and x by x llk , we obtain 

_ x 

y " 1 — kai, k +ix' 

Replacing next, y by \jy, and x by 1/x, we have 

y = x — ka l)k +\, 

so that in this way we have secured a transform of y = A x (x), in the form 

y = c + x. 
Putting — kai, k+iy for y and — Mi, k+iX for x, we arrive at exactly 

y = 1 + x, 

as desired. We need merely retrace our steps to secure the transforming 
series B(x). 

In the case of Type He, the best that we can do for the first two steps 
is to transform y = A\(x) into 



V 



1 — kai, k+ix — k ai i2 k+i ?> — ^ k + 1 '• 



This is an analytic function in the neighborhood of the origin, that resembles 
somewhat the analytic function 

x 
y ~ 1 - ex' 

where c is a suitably chosen constant. Any transforming function, however, 
which carries the given function into one of the latter type, must have an 
essential singularity at the origin, it being understood that for Type lie, 
25i, 2k+i 4= (Jo + l)5i 2 , jb+i. For a fairly complete discussion of the existence 
of this transforming function, cf. Leau, loc. cit. It will be unnecessary for 
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us to go into this particular question, since a direct representation in series 
may be found without the use of a transforming function. We need only 
regard Case II, as a limiting form of Case I, and pass to the limit in either 
(4) or (5). We obtain in both cases the same series.* We shall have: 

U (x) = A (x) = x, 

Ui(x) = Ai(x) - A (x), 

U t (x) m A 2 (x) - 2A 1 (x) + A (x), 



Ui{x) s Ut-AAiix)] - ?7_i(x), 

= At{x) - \l J Ai-fa) + (ijAi-iix) - ••• ± A (x), 



while for on we have 

n — i + 



-«,-! = (g), 



Hence, when a u — 1, ^- ^., for Case II, 

(6) A n (x) = A + jUi - A ) + n{n i 7 2 1) ^ - 2A > + A «) + • • •• 
For Case III, if Ai(x) be given, there exist series C(x), for which 

C[A 1 (x)] m A.ldx)], 

but these C's, instead of containing but a single parameter, contain more 
than one, and in some cases, contain an infinite number of parameters. 
For example Ai(x) = — x is commutative with any odd function of x. 
But if Ci(x) be commutative with A\{x) then C m (x) is commutative with 
A n (x), where m and n are any positive integers, although the converse is 
not necessarily true. For c n an integral root of unity, we may choose m 
and n so that a m \ = 1, and c OTi = 1. This carries us back to Case II. In 
particular, by taking m = 1, we may conclude that if A\(x) be such that 
a Ci(x) exists for which c n is equal to one, then Ci(x) can contain, while 
cii is kept equal to unity, at most one free parameter, unless Ai(x) is of 
Type Ilia. Those of Type Ilia, all contain an infinite number of parameters 
as we shall presently show. The same remarks apply when Cu is any 
integral power of a n , say the Tth, since Ci[A h -i(x)] will also be commutative 
with Ai(x), where h is the lowest positive integer for which n n h — 1, but 
Ci[A h -i(x)] will have unity for the coefficient of the first term. Apart 

* This series has been obtained directly by Schroeder, Math. Ann., 2 (1870), p. 317. 
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from Cu, the only arbitrary parameter is a certain Ci,*+i, as is readily 
verified by the method of undetermined coefficients. Now x itself, is 
certainly commutative with Ai(x), so that there exists at least one series 
Ci(x) for which c n = 1, cu = 0, i = 2, 3, • • •. Since all series Ci(x) with 
the same Cu coincide as far as the term containing the parameter C\, k +i, 
it is clear that A h (x) which is commutative with A\(x), and for which 
ahi = 1, is of the form 

A h (x) = x + 0+---+0 + a h , k+1 x k+1 + a h , k+2 x k+2 + ■■■. 

The method of undetermined coefficients which gives the k to be used, 
shows us that this k is h itself, or an integral power of h. Hence had we 
started with a series 

A 1 (x) =x + 0+---+0 + ai, k+1 x k+1 + oi, k+2 x k+2 + ■■■ 

where ai,jt + i 4= 0, then the only expressions A a / h )(x), which can exist, will 
be for h, a factor of k. For types III& and IIIc, integer iterates, of course, 
always exist, and they will either be of the same type as the given series 
or of the analogous type under Case II. In some instances a few fractional 
iterates of the type of the given series may exist. However, no series A n (x) 
of Types Ilia, or Illb, can exist, containing an arbitrary parameter, n. 

Type Ilia demands further attention.* For the sake of definiteness, 
we shall first consider the case in which h = 2. Here A x (x) is such that if 
y = Ai(x) then x = Ai(y). Viewed geometrically, this gives us a curve in 
the x, y plane, whenever the series A x (x) is convergent. This curve is 
unaltered by interchanging x and y, that is, it is symmetrical about the line 
y = x. If we transform the coordinates by writing 

u Q = JOc - y), 

Mi = \{x + y), 

the curve will be symmetrical about the line u = 0, i. e., the Mi-axis, 
and in the neighborhood of the origin it will have contact with the M -axis, 
since the slope of the curve y = Ai(x) is unity, at the origin. Thus, if in 
the neighborhood of the origin, u x be represented as a power series in u , 
for points on the curve, the series will represent an even function of u . 
Furthermore, this is the only restriction imposed upon u u as a function of u . 
We may select for U\, an arbitrary even function of u , and determine 
x and y by the equations 

X = U + Ui, 

____^^_ y = — u + Mi. 

* This leads to Ah(x) = x, which is known as " Babbage's Problem." Cf. the following : 
J. D. Gergonne, Ann. Math., 12 (1821-2), p. 73. O. Rausenberger, Math. Ann., 18 (1881), 
p. 379. E. Iaggi, Nouv. Ann. (3), 19 (1900), p. 483. M. J. van Ufen, Amst. Ak. Versl., 
18 (1909), p. 860, and 19 (1910), p. 27-31. The discussion here is new. 
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If from these equations we solve for y as a function of x by eliminating the 
u and Mi, the solution y = A\(x), will be such that also x = A\{y). We 
may forego the geometric picture, and select Ui as an arbitrary series in 
even powers of u without constant term, and not necessarily convergent. 
The series A x (x) obtained by solving as before, for y, will be determinate, 
and if y = Ai(x), x = A\{y). 

The discussion for a general integer value of h, follows along similar lines, 
although no simple geometric illustration of an obvious sort suggests itself. 
We shall first suppose a series A\(x) given for which Ah(x) = x. From this 
series we shall determine h series Ui, i = 0, 1, • • •, h — 1, which we shall 
find to be restricted in a very simple manner only. This done, we shall 
start with the h series u t and so determine the form of the most general 
series Ai(x) for which A% = x. Let us write 

Wo = ^ (A + wii + oM 2 + • • • + o^Ak-i), 
U! = I (A„ + w 2 Ai + W <A S +■■■+. w 2A - 2 A A _i), 



w/,-1 = ^ (A + Ai + A 2 + ■■■ + A h -!), 

where wis a primitive hth. root of unity. We wish now to discover what 
restrictions must obtain among the u's in order that Ah = A . If we 
replace A,- by A,-+i, Mo is carried into w'« , where w' is the reciprocal of w, 
while ttA_i is unaltered. Hence if we express Uh-\ as a power series in u , 
by eliminating the independent variable x, occurring in each A,-, i = 0, 
1, • • •, h — 1, we shall find that only the terms in powers of u h appear, 
for if any other term were to appear, it would be changed into the product 
of its original value by a power of w, different from unity, which is impossible, 
since u h -i is to remain unaltered. In a similar manner the power series 
for Ui, i = 1, 2, • • •, h — 1, in terms of u , contain only terms in which 
u i+1 times a power of u h appears. These are the only restrictions upon the 
series «,- in terms of u . We may suppose, if we desire, that the u's are 
given, and that we are to find Ai(x). We may take any h — 1 power 
series, Pi{t), i = 1, 2, •■•, h — 1, convergent or divergent and write 
Ui = Uo i+1 Pi(u h ). If we solve for the A's formally in terms of the u's, 
we obtain in particular 

A (X) = X = (tt + Ml + • • • + Mft-l) 

i. e., a; = E-i(uo) where E-i(u ) is, by definition, the series in w obtained 
by replacing each Ui, i = 0, 1, 2, • • •, /i — 1 by its series representation in 



THE ITERATION OF FUNCTIONS OF ONE VARIABLE. 39 

terms of u , in the right hand member of the above equation. By definition 
of E-i, it follows that 

1*0 = E(x). 
Now E_i(m ) is of the form 

JE_i(tto) = wo + i*o 2 Pi(i*o A ) + M 'Pi(Mo*) + • • • + i*o A P*-i(i*o*). 

Thus E_i(o)i* ) is simply (wi* + w 2 i*i + w 3 i* 2 + • • • + i*&-i), since when 
we replace 1*0 in it* = i*o t+1 Pt(i*o A ), by wuq, we obtain merely co i+1 i* t -. But this 
expression for E-i(a)Uo) is exactly what we obtain formally, when we solve 
for A\ in terms of the i*'s. Thus we have A\ s= E-i(cou ) or 

A x (x) = ^[co^x)], 

from which we verify that indeed A h (x) = x, since obviously for i any 
positive integer we have 

Ai(x) = B_i[«*^(x)]. 

Thus we verify that it is not only necessary but sufficient for i*;(i* ) to be 
of the form Uo i+1 Pi(u h ), in order that the Ai(x) associated with the i*'s 
as above shall be such that A h (x) = A (x). The E-i(u ) that we have 
obtained above is a power series in i* characterized by the property that 
the absolute term is missing, the coefficient of u is unity, and the coefficient 
of u lk+1 , I = 1, 2, • • •, vanishes. For any given series Ai(x) the series 
E-i satisfying the conditions just mentioned is uniquely determined by the 
equation 

Ai(x) = E-i[aE(x)] 

analogous to the equations obtained in Case I. Unlike the situation in 
Case I, we have here imposed additional conditions on E-i. If we replace 
1*0 by a power series 1*0 = S-i(v), where S-i is restricted to the extent 
of having the coefficient of every term which does not contain v to a power 
of the form Ik + 1, I = 0, 1, 2, • • •, vanish, while the remaining coeffi- 
cients are arbitrary, then we obtain 

A<(x) = E^iS^SiEix)}]}, 

where A it E, and o> are fixed and *S_i has an infinite number of arbitrary 
coefficients as already remarked. We may think of E-i[S-i(v)] as con- 
stituting a new series P_i(i>), for which the coefficients of v lk+1 , I = 1, 2, 
are arbitrary, the others being then determined, if A i be given, so that 

Ai(x) = F.i[aF(x)]. 

Case IV presents little of interest. Type TVa is trivial, and will be 
passed over without further comment. Type IV6 has a general iterate as 
a function of n, but it is not expressible as a power series in integral powers 
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of x; terms with irrational exponents, for example, must be introduced. 
We may however obtain a transforming function, which reduces this case 
to Case I. We first choose a, as one of the (k — l)st roots of au, and 
replace y by y/a, and x by x/a. We obtain in this way a new series of 
Type TVb, but one in which a ik is now unity. By the method of undeter- 
mined coefficients, we may find a series F(x) such that 

F[Ax(x)] - [F(x)]\ 

so that Ai(x) is now transformed into x k . The equation y = x k is trans- 
formed into y = kx, by replacing y by e v and x by e x , while y = kx is of the 
form Type 16. 

The Transforming Series as a Limit Series. 

8. We shall now suppose that an Ai(x) and a Gi(x) are given and known, 
while the existence of an Fi(x) is known but not the form of the series, where 

Ai(x) mFiGiiF.^x)]). 

We shall also suppose that the general iterate G„(x) is known, and we 
wish to determine F so that 

A n (x) = F.{Gn[F-i(x)]}. 

We shall now describe a method which in some important cases serves to 
define an F{x). The method is not always applicable but is of considerable 
interest in the cases in which it can be used.* 

For any given positive integer p, let us consider the four series G- p [A p (x)], 
Gp[A- p (x)], A- p [Gp(x)], A P [G- P (x)], of which the latter two are the in- 
verses respectively of the former two. A given term in any one of these 
series will be a function of p, which may or may not remain finite as p 
takes on successively larger and larger positive integer values. If it should 
happen that in one of these series the coefficient of each term approaches 
a definite finite limit for p = + a, then a limit series, L(x) will be defined. 
The question is one of convergence, not of the series as a series in x, but of 
the coefficients of the series as functions of p. If L(x) exists, its inverse also 
exists, so that we may without loss of generality suppose L to be the limit 
series of one of the first two mentioned above. From the definition of 
L(x) it then follows that 

G[L(x)} = L[A(x)\, 

so that L(x) is a series of the form F-i(x) which was to be found. In 
determining L{x) we have not made use of the actual form of G n {x), for n 
other than an integer, so that when different methods of interpolation are 
used to give different forms for (?„(x), correspondingly different forms will 

* Cf. G. Koenigs, Ann. Ec. Norm. (3), 1 (1884), Suppl., p. 19, and (3) 2 (1885), p. 385. 
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be obtained for A n (x), n other than an integer. If G- p [A p (x)) and G P [A- P (x)] 
both have limit series L, and these are distinct, then for every integer value 

of n, we shall have 

G n [L(x)] s L[A n (x)] 

true, for both choices of L(x), but for fractional, irrational or complex 
values of n, it may well happen that the same determination of G n (x) 
gives for the two L's different determinations of A n {x), which corresponds to 
the arbitrary feature of G n (x). If Cix) be any series commutative with 
G{x), then C[L(x)] may be used equally well, in place of L{x). In particular 
C(x) may be an integral or any fractional iterate of Gix). But this corre- 
sponds merely to having considered G- p+q [A p (x)] or G p - q [A- p (x)], respec- 
tively, in place of the given series, where q is a fixed number. Similarly 
we might have taken L[C(x)] in place of L(x), where C{x) is this time com- 
mutative with A{x). These two cases are not, of course, independent, 
since, as we have already remarked, G„[L{x)] = L[A„(x)]. The one- 
parameter family of series G g [L(x)], q arbitrary, which exists if L(x) does, 
constitutes in many cases the totality of series F_i(x). 

Investigation of the Convergence of the Series Employed. 

9. We now come to the question of convergence. In Type la, whenever 
A\(x) is a convergent series, the convergence of E(x) follows, and con- 
versely, as may be proved in several ways.* Type 16, although not usually 
treated, presents no difficulties, since Ai(x) and A_i(x) yield the same E(x), 
and since when \a n \ > 1, then | a_i,i | < 1, a-i,i being merely 1/an. 
Thus, if Ai(x) is convergent and | an | > 1, then E(x) is convergent, 
because of the convergence of A-i(x). Type Ic is of a different sort. For 
A\(x) convergent and of Type Ic, E(x) may be either convergent or diver- 
gent. If Ai(x) is convergent certain conditions must be satisfied in order 
that E(x) be convergent. These are not formal conditions upon the first p 
coefficients, p = 1, 2, • • •, since the problem is one having significance only 
for convergent series ^i(x), and convergence is not dependent upon the 
values of any finite number of coefficients. If E(x) be convergent, A(x) 
must be. Necessary and sufficient conditions for the convergence of Eix), 
Ai{x) being given as convergent, may be given in different forms. We 
shall state a condition which theoretically is capable of explicit application 
in any case. The cases of Type Ic, for which E{x) is convergent, we shall 
call of Type Ic', those for which E{x) is divergent, being of Type Ic", 
provided that Ai(x) is convergent, otherwise of Type Ic'". We shall now 
state 



* For example, Leau treats the equation in this instance both by the method of dominant 
functions, and by that of successive substitutions, loc. cit. 
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A necessary and sufficient condition for Type Ic'. It must be possible 
to find a circle about the origin in the complex plane such that if any point 
P other than the origin be taken within this circle, the totality of the 
positive and negative iterates of P, constitute a set of points lying every- 
where densely upon a simple closed curve of finite length, and with con- 
tinuously turning tangent, encircling the origin. 

This includes, in particular, the conditions of interior and exterior 
"stability," about the origin. The condition of interior stability, frequently 
called merely "stability," is the following: For any sufficiently small circle 
Ci with center at the origin, there must be a smaller concentric circle Ci, 
such that the positive and negative iterates of every point within C 2 is 
within Ci* The condition for exterior stability insures that every iterate 
of a point in the neighborhood of the origin but exterior to a small circle 
enclosing the origin remains exterior to the same or a smaller fixed circle 
about the origin. 

To prove that the condition stated above is sufficient, we may notice, 
that if it be satisfied for a given P, we shall have a curve C completely 
determined by the iterates of P, together with the conditions that C is 
continuous, and its tangent turns continuously. This curve C may be 
transformed into a circle in the plane of a complex variable t by a transforma- 
tion of the form t — Z?_i(x), where E-i(x) is an analytic function, as follows 
from the theory of conformal mapping. The function E- X will be com- 
pletely determined when we require the origin in the x-plane to go into the 
origin in the £-plane, the circle in the i-plane being the unit circle and the 
positive real axis in the x-plane going into some sort of curve tangent to 
the positive real axis in the t-pl&tte, at the origin. Any analytic trans- 
formation of the x-plane in the neighborhood of the origin, which carries 
the curve C into itself and the origin into itself, will correspond analytically 
to a transformation of the 2-plane in the neighborhood of the origin, leaving 
the unit circle and the origin unaltered. Now Ai(x) is such a transforma- 
tion of the x-plane, and it corresponds therefore to a rotation in the i-plane, 
since a rotation is the only conformal transformation of the unit circle into 
itself which leaves the center invariant.! But formally it corresponds to a 
rotation an-t, | an | = 1, where a u is the initial coefficient of the series 
Ai(x). Thus the formal series E{x) is actually the inverse of the series 
corresponding to the analytic function 2£_i(x) just determined. Hence the 
above condition is sufficient; that it is necessary is seen immediately from 
the fact that if E is analytic then the successive iterates of a point P in the 
x-plane correspond to the successive images of a point in the £-plane, as the 
2-plane is rotated through an angle incommensurable with 2ir. 

* Cf. Levi Civita, Ann. di. Mat. (3), 5 (1901), p. 240, and Cigala, Ann. di. Mat. (3), 11, 
(1905), p. 67. . 

f Cf. Osgood, Funktionentheorie, vol. I, 1st ed., p. 595-. 
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In the case of Type Ila, we may transform x into itself by the use of 
divergent as well as of convergent series. No simple or direct proof of the 
convergence of the series transforming Type lib into 

x 

y 



1 — 2a,i 2 x ' 



seems to have been made. The convergence follows, however, when Ai(x) 
is itself, convergent, from the detailed investigation of a function-theoretic 
character made by Leau.* Since the series used in Type He are not 
formally equivalent to 

x 
y ~ T^W^z' 

no convergent power series can be found which transforms the former into 
the latter. The existence of functions in general analytic, but possessing 
essential singularities at the origin is proved by Leau for the cases in which 
the given series is convergent. 

The existence of divergent series of the Type Ilia has already been 
shown. Whenever the series is convergent, the h series, Ui, i = 0, 1, • • • , 
h — 1, are convergent, since they are linear combinations of the first h 
iterates of Ai(x), so that the corresponding E(x) must be necessarily con- 
vergent, while for A i(x) , divergent, E(x) is divergent. In the cases of Types 
Ilia and Illb, for which, in general, no iterate containing an arbitrary 
parameter n, can exist of Case III, it is not surprising that no transforma- 
tion is known which reduces either of these types to the form ex, or x + a , 
and certainly no conformal transformation of this sort can exist, so that 
there is no question of convergence. 

A discussion of the convergence of the transforming function E, which 
reduces Type IVb, to the form y — x k has been made by Grevy, loc. cit. 

A Geometric Interpretation of Iteration. 

10. We have already mentioned an infinite matrix whose first principal 
minor of order to, we have denoted by Ai (m \ which may be considered as 
serving to define the series Ai(x). We have mentioned the problem of 
the reduction of this infinite matrix to a normal form by successively 
reducing to normal form the matrices Ai (m) , to = 1, 2, • • •. In place of 
the collineation symbolized by y — Ai (m) x, and its reduction to a normal 
form, we may consider the entirely equivalent problem of the simultaneous 
linear transformation of the pair of bilinear forms, or of the X-matrix, 
corresponding to the original matrix. In terms of the original series, Ai(x), 
this means that we may consider the analytic transformation of the pair of 

* Leau, loc. cit. 
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expressions y = x, and y = A\(x), or of the linear function of X, expressed 
by y = A\(x) — \x. Since any analytic curve in the real plane may be 
transformed analytically in the neighborhood of a non-singular point of the 
curve into y = x, the invariants under analytic transformations of a pair 
of curves at a point of intersection, are the invariants of the pair y = x, 
and y = Ai(x), where A\(x) is the transform of the second curve. This 
excludes, of course', certain exceptional cases, as when a derivative of Ai(x) 
becomes infinite, which, however, may be readily treated as limiting cases. 
The group of transformations, y' = <p(y), x' = <p(x), has significance 
chiefly in the " function-plane," and by change of variables, in the "inver- 
sion-plane." The study of the conformal invariants of a pair of real curves 
at a real point of intersection, is essentially the study of Cases II and III, 
discussed above, while the general conformal group without regard to 
reality involves also Cases I and IV.* Integral iteration corresponds to 
the multiplication of the angle between y = Ai(x) and y = x, while division 
of the angle is secured by fractional iteration. 

PART TWO. THE ITERATION OF A REAL FUNCTION. 
The Domain of Definition of A n . 

11. We shall consider the iteration of a real function defined for the 
points within an interval a, b, a < b, where, in particular, a may be nega- 
tively infinite, or b positively infinite, or both. For many problems in 
analysis, it is immaterial whether we consider a real analytic function as 
defined also for imaginary points, or not, so far as the phenomena for real 
points alone are concerned. Iteration is not, however, one of these prob- 
lems, since it frequently happens that the iterate of an imaginary point is 
itself a real point. 

Let Ai(x) be a finite, real, one- valued function of x, defined over a range 
ri in the interval a, b. We shall usually suppose that the range r/ of values 
taken on by x' = A\{x) lies within the interval a, b, since otherwise, we 
might have extended our original interval a, b, so as to include r/, without 
altering r x or any other data. 

The function A 2 (x) = AilAtix)], will be defined only for the numbers of 
the subset r 2 common to r x and r/. If we denote by r 2 ', the values assumed 
by A\(x) as x varies over r 2 , then A 3 (x) is defined only for the subset r 3 , 
common to r 2 and r 2 ', and so forth. It may happen that after iterating a 
finite number of times, we may eventually find an n, so large that r n is a 
null-set, and no expression A n (x) can be defined. This case, although 
common, presents but little of interest. If n be an open set, it may happen 

* Cf. Kasner, "Conformal Geometry," loc. cit., and Pfeiffer, notice in Bull. Am. Math. 
Soc, 21 (1915). 
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that for every positive integer n, r n contains some number, while no number 
is common to all r n 's, owing to the fact that the limits approached lie out- 
side of ri. A third case arises when the limit set, r a , consists of a finite 
number of points, e\, e 2 , • • • , e&, so that these are the only points for which 
A n (x) can be defined for every positive integer n. Now Ai(e t ) must be equal 
to one of the e's, say ej, for if A n (ei) is denned for all positive integral values 
of n, so must also A„_i[.Ai(e;)], since these are identical by definition, or in 
other words -Ai(e,), is an e by definition. Since Ai(x) is one-valued, each 
e» must go into a unique e ; - either itself or another. Of the k, e's, e h e 2 , 
■ ■ • , ek, it may happen that several distinct e's are carried into a single one 
by the operation x' = Ai(x). There must be a non- vanishing subset of 
these k, e's into each one of which subset, an e of the original set is carried 
by x' = Ai(x). By repeated selections of such subsets, we obtain finally a 
non-vanishing set of k', e's, which we may call e\, e 2 , • • • , eu>, such that each 
e of this set is carried into an e of this set and is obtained from an e of this 
set by any transformation x' = A n (x), n, a positive integer. Since k' is 
finite, and each e, is carried into but a single e,, each ej must be obtained 
from but a single e, of the set of k' , e's. We have then a permutation group 
upon k' elements, and the transformation x' — A t > t (x) must leave each of 
this set unaltered. Other cases arise in which r K coincides with r\. It may 
possess no fixed points, a finite number, or an infinite number, and every 
point of r K may be fixed. 

We may distinguish in all cases, as we have already done in one case, 
a subset of r x , which we shall call rj, where not only is A n {x) defined for 
every point of rj, for n an arbitrary positive integer, but A-„(x) is also 
defined for all positive integral values of n. For the subset of r x obtained 
by removing rj, we have points which are carried into other points by A n , 
but which are not obtainable from other points by an A n , where n is a 
sufficiently large positive integer, in each case. The subset denotable by 
?"» — rj, may be further subclassified. Regarded as a set of points it is 
completely arbitrary in character. The set rj, may be divided into sub- 
sets, r xl ', r^, • • • , together with certain residual sets r a \ eU , rj p0 , r and r K \ hyp . 
The set rj p , for p a positive integer, is defined as a set which is carried 
pointwise into itself by A P (x), but by no A„(x), n < p. The set r„\ M is 
carried into itself circularly, but no point is carried into itself exactly by 
any A n , n, a positive integer. This corresponds to the case of a rotation 
through an irrational angle. The set rj t hyp , is a set carried into itself by 
A n (x), n a positive integer, but not pointwise, there being for any par- 
ticular case two limit points, after the manner of a hyperbolic linear 
fractional transformation. Similarly any particular set rJ iPax has one 
limit point. The set r Kl ' is absolutely arbitrary in form. The set r wl ' is 
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restricted to the extent that it must consist of a set of pairs of points, 
without repetition, or as it may be regarded, as a pair of sets of points, 
the points in one set being in one to one reciprocal correspondence with 
those of the other. Similarly the set r xn must consist of sets of n points 
without repetition, or n sets in one to one reciprocal correspondence with 
one another. The sets, r a \ ell , may be regarded as constituting independent 
systems, as also in the cases of rJ iPIIS and r x \ ell . 

Some General Theorems on A n . 

12. We shall state a few of the numerous theorems which are valid 
for real one-valued functions defined over a range r in an interval a, b, a < b, 
extending, in particular cases, to infinity. The method of proof is fairly 
obvious in every case, and as the proof itself is never difficult, it will be 
omitted. The index n, of iteration will be restricted to positive integers. 
We shall use the notation c, di for an interval, c ^ x ^ di included within 
the given interval. 

If Ai(x) is periodic, A n (x) is periodic in so far as it is defined. 

Every point of intersection ef y = Ai(x) with y = x, is also a point of 
intersection of y = A n (x) with y = x. 

If Ci ^ Ai(x) ^ di, whenever c 2 ^ x ^ d 2 , for c 2 ^ Ci < di ^ d 2 , then 
the same is true for A n (x). 

If :s si < s 2 , and Ai(e) = e, and if also for x in the interval c, d, where 
c < e < d, A x (x) lies between e + Si(x — e) and e + s 2 (x — e), then for x 
in the same interval, A n (x) will lie between e + Si n (x — e) and e + s 2 n (x — e) 

If e be a value for which Ai(e) = e, then every point of intersection of 
y = Ai(x) and y = e is also a point of intersection of y = A n (x) and y = e. 

Hence we obtain also the following theorem, 

If c ^ Ai(x) ^ d, for a ^ x ^ 6, and if y = x for c ^ x ^ d, then 
A n (x) = Ai(x). 

If for every Xi < x 2 , Ai(x) < ^1(2)2, then also A n (xi) < A n (x 2 ), in so 
far as A n (x) is defined. 

If the absolute minimum of Ai(x) is unique and occurs at c, and Ai(c) > c 
and also the absolute maximum of A i(x) is unique and occurs at d, and 
Ai(d) < d, and if Ai(x) in so far as it is defined is an increasing function for 
c ^ x ^ d, then whenever A n (c) is defined it will be the unique absolute 
minimum value of A n (x), and A n (c) > A n _ t (c), i = 1, 2, • • •, n, in so far 
as these are defined. A corresponding situation holds for A n (d), and 
A n (x) will be, in so far as it is defined, an increasing function of x for 
c ^ x ^ d. 

If for every Xi < x 2 , 4i(xi) > Ai(x 2 ), then also A 2n -i(xi) > A 2n -i(x 2 ), 
while A 2n (,xi) < A 2n (x 2 ) in so far as A 2n ^i(x) and A 2n (x) are defined. 
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If A\(d) is either undefined or is equal to Ci and for a == Xi < x 2 — c 2 , 
Ci ^ Ai(xi) < Ai(x 2 ) == 6, and if for c 2 ^ £1 < £2 — &, a ^ Ai(xi) < 
Ai(x 2 ) ^ Ci, while Ai(b) is either undefined or equal to or less than Ci, 
then there exists a pair of number Ci (n) and C2 (n) , for which an analogous 
situation holds for A n (x). 

If in the preceding theorem we put "s" in place of "<", and 2n — 1 
for n, the theorem continues to hold, but if we put "s" in the hypothesis 
and 2n in the conclusion, we must read "s" in the conclusion. 

Other theorems might be mentioned of a nature analogous to those just 
cited. We shall now however give a few of a different type, in which the 
continuity of Ai(x) is supposed, as well as the existence and in some cases 
the continuity of the first derivative of Ai(x) with respect to x. The proofs 
are obvious and will be omitted. They involve chiefly the use of the 
identity, 

A n '(x) = A 1 '[A^(x)]-A l '[A n -i(x)]---A 1 '[A l (x)]-A 1 '(x), 



where A n '(x) denotes 



9A n {x) 

dx 



If the total number of maxima of Ai(x) be finite and equal to m, while 
Ax'(x) does not remain zero over any interval, the least possible number of 
maxima of A n (x) will be m — 1, and the greatest possible number will be 
2 n ~ 1 m n . Examples of the two cases are most readily grasped by means of 
figures. Fig. 1 represents a broken line joining in order the points whose 




Fig. 1. 



coordinates are, [0, l/(m — 1)], [l/(m — 1), 0], [2/(m 
[3/(ro - 1), 0], • • -, [2 - l/(m - 1), 0], [2, l/(m - 1)]. 



1), l/(m - 1)], 
Odd iterates of 




Fig. 2. 
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this curve coincide with the curve itself while even iterates are of the form 
illustrated in Fig. 2, where we have a broken line joining in order the points 
whose coordinates are [0, 0], [l/(m — 1), l/(w — 1)], [2/(m — 1), 0], 
[3/(m-l), l/(m-l)], .-., [2-l/(m-l), l/(m - 1)], [2, 0]. An 
example of the second case is given in Fig. 3, by the broken line joining in 




A/(x) 

throughout the interval, 
tion, so also is y = A n (x). 



Fig. 3. 

order, the points, [0, 0], [1/ro, 2], [2/to, 0], [3/m, 2], [4/m, 0], • • •, [2 - 1/m, 2], 
[2, 0]. The above theorem holds also if " minimum " be read for " maxi- 
mum." 

If for every x, \Ai(x) \ < s, where s is a constant, than | A n '(x) | < s n , 

Hence if y = Ai{x), be a curve of limited varia- 
The length of arc of y = A„(x) between x = c 
and x = d, c < d, will be less than (d — c)s n . 

If < Si < Aiix) < s^ then < Si n < A n '(s) < s 2 n . 

If | Ai'(x) | < 1, throughout the interval, then the number of maxima 
of A n (x) cannot differ by more than one from the number of maxima of 
Ai(x), provided that Ai(x) does not remain equal to zero over an interval. 
The same theorem holds true if " minima," be read for " maxima." 

If throughout the interval, Ai(x) ^ a + (x — a) k /(b — a) k ~ x , then 
A n (x) < a + (x - a) kn /(b - a) k "~ 1 . The same theorem holds if "<," be 
replaced by ">." 

One-Valued Iteration, A n , for n Varying Continuously. 

13. Let us examine those cases in which Ai(x) is one-valued, and A n (x) 
exists as n varies continuously over real values. We shall have given a 
range r over which Ai(x) is defined. Inr, we shall distinguish two classes 
of points, which we shall call singular points, and ordinary points. An 
ordinary point will be defined as a point of r, for which there is no next 
point of r whether to the right or to the left of it, the other points are the 
singular points. A point which is not of r, but which is a limit point of 
points of r, and for which on both sides, there is no next point of r, will be 
called a gap point of r. We shall suppose that except in the neighborhood 
of x — a, and x = b, A n (x) exists for all values of n for which | n \ is suffi- 
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ciently small. Since for n irrational, A n {x), is defined only by means of 
interpolation, we shall be forced to regard A n (x) for any given x in r, as 
varying continuously over r, when n varies continuously over real values. 
This does not mean that A n (x) for a given x varies continuously in the 
interval a, b, as n varies continuously, but simply that it varies continuously 
with reference to the inner structure of r. Since continuous variation from 
a singular point is impossible, we must have A n (x) = x, for every singular 
point. It must be borne in mind that a point is singular or not, purely in 
terms of the internal structure of r, and the end point of an interval c, d of 
points of r need not be a singular point. Since by continuous variation in r, 
a gap point may be approached but not passed, we see that the limit points 
approached by A n (x), for x given, and n increasing by integers toward 
infinity include in general gap points. If we take three numbers, x x < x 2 
< x 3 , then for \n\ small, we must have A n (xi) < A„(x 2 ) < A n (x s ), pro- 
vided, at least, that a < xi and x 3 < b, since a continuous group of trans- 
formations of the form here considered cannot include a transformation 
by which sense is altered, except possibly in the neighborhoods of a and b, 
for | n | small. We may make three possible conventions with regard to 
the end points a and b, which will be in accord with the conditions imposed 
upon interior points. 

First. We may require A n (x) to be defined at a and b only when n is 
positive for one of these points and negative for the other. 

Second. We may require A n (x) to be defined at both a and b for [ n \ 
small, and for n both positive and negative. 

Third. We may regard the interval a, b, as treated cyclically, so that 
a and b are not distinguished, and a point to the right of a is treated as to 
the right of any point to the left of b, only cyclic order being essential. 
In this case, we may, if desired, require, in particular, that r coincides with 
the whole interval a, b, and A n (x) is defined for every point of the interval, 
whenever | n \ be small, whether n be positive or negative. 

In the first case, a and b will be essentially different from any other 
points of r, if they, or either one, be in r. The range of definition for A n (x), 
will in general depend upon n, and various small difficulties arise. It is 
however a common and important case. In the second case, a and 6, if in 
r, are merely singular points, as already defined, and are not to be dis- 
tinguished from the other singular points. In the third case, a and b may 
be regarded as a single point in no way distinguishable from any ordinary 
point or r, either in terms of the structure of r itself, or in terms of the 
behavior of A n (x), as a function of n. 

We may describe the curve y = Ai(x), where Ai(x) satisfies the condi- 
tions of this section as being with respect to r, a series of monotonically 
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increasing curvilinear segments, and discrete points. The discrete points 
are all on the line y = x, and have for their abscissas, the singular points 
of r, measured on the x-axis. The curvilinear segments, if they have end 
points, have these on the line y = x. If we require that the range of defi- 
nition interior to a, b, with the possible exception of the neighborhoods of 
a and b themselves, shall not vary with n, then we have also that curvi- 
linear segments, where there is no end point, approach, i. e., at a gap point, 
the line y = x, and have a point of y = x for limit point. When referred 
to the interval a, b, without the intervention of r, a more complicated 
description of y = A x (x), is, of course, necessary. 

The Transforming Function B(x). 

14. In defining A n (x) for n, a fraction, we may make use of a trans- 
forming function, B(x), such that if y = A„(x) 

jx = B(t), 

\y = B(n + t), 

so that A n (x) = B[n + B-i{x)\. If y = A n (x) is to be a one-valued func- 
tion of x, for positive values of n, then B(t) must be monotonic if it be con- 
tinuous, that is, if A„(x) be a continuous function of n for a given x. For 
suppose if possible that B{t) is one- valued and continuous but not mono- 
tonic. Then let B(a + h) = B{h), whence since A„[B(t x )] = B(a + h), 
we have A a [B{ti)} = B(h) and therefore A ma [B(ti)] = B(ti), where m is 
any positive integer, or B[ma + ti] = B(U), m = 1, 2, • • •. B(t), being 
supposed non-monotonic, will have either a maximum or a minimum, if not 
many of each. Let us suppose there is a maximum /3 = B(y). A line 
x = /3 — «, e > 0, will cut the curve, x — B(t), at two points in the neigh- 
borhood of (7, 0), which we may call t = y — Si, and t = y + S 2 , 5i > 0, 
§2 > 0. Making use of our previous remark, we conclude that 

B[m(h + 81) + y - 5J = B(y - 50, m = 1, 2, • • •. 

Thus the line x = (3 — e contains a set of points of the curve x = B(t) 
beginning with t = y — 8 h and succeeding each other at a constant interval 
Si + 5 2 . As e is taken smaller and smaller, we obtain new lines and each 
contains points of the curve x = B(t) distributed in the same fashion at 
constant intervals from each other, the constant interval, however, decreas- 
ing with e. From the continuity of B we would conclude that in the limit- 
ing case of the line, x = /3 every point of this line to the right of t = y, ought 
to be a point of the curve x = B(t). But if x = B(t) is to be one- valued, 
such a situation is impossible. We therefore conclude, that for x = B(t), 
continuous and one- valued, this curve is also monotonic. The proof for 
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the case of a minimum but no maximum follows in the same fashion and 
is therefore supposed completed. 

The range of the variable t, has, of course, nothing to do with the range 
r or the interval a, b. If B(t) is mono tonic and one- valued, its most sig- 
nificant characteristics will be its vertical and horizontal asymptotes if 
one or both of these exist. If u = B(t) has a horizontal asymptote u = e, 
then x = e, is a fixed point for y = A n (x). If u = B(t) has a vertical 
asymptote, then y = A n (x) has for every n, a fixed common horizontal 
asymptote, and the interval a, b must be infinite, while if the curve u = B(t) 
be defined to the left of a vertical asymptote, then y = A n (x) has a vertical 
asymptote. Further theorems of a like sort are readily obtained, for in- 
stance: If u = B{t) approaches minus infinity, as t approaches e, from the 
left, then y = A n (x), is an increasing function with a vertical asymptote 
x = e — 1, and lying to the right of the asymptotes. 

A transforming function Bit) must be itself monotonic and is adapted 
for the study of only monotonic functions y = Ai(x). If in the neighbor- 
hood of a point x = y = e, y = Ai(x) is monotonic and increasing, we may 
obtain for this neighborhood a transforming function x = B{t). The 
question arises as to whether if we continue Aiix) beyond the neighborhood 
in which it is monotonic, there will not be a corresponding extension of 
Bit), so that the B(t) so extended shall still serve as the transforming func- 
tion of y = Ai(x) when this latter is extended. We shall find that no ex- 
tension of B(t) can be made to serve as a transforming function for Ai(x) 
when extended beyond the region in which A\{x) is monotonic. If, how- 
ever, Ai(x) is analytic, we may determine B(t) analytically, and each will be 
continued analytically beyond the region in which it is monotone. What 
then is the relation between these continuations? For y = Ai(x), one- 
valued, but not monotonic, the extension beyond the monotonic interval 
containing a point x = y = e, leads us to an extension of B(t) into an 
auxiliary part, such that x and y may both vary on the principal part B(t) 
with the relation x = B(t),y = B(l + t), but x alone varies on the auxiliary 
part x = B(i), while y continues to satisfy y = 5(1 + 0- The situation 
may readily be grasped upon inspecting the case in which A(x) = sinx. 
The curve u = Bit) is indicated, Fig. 4, by a heavy curve, while the 
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auxiliary part u = B(t), here consists of an infinite sequence of branches 
a few of which are here indicated by dotted curves. 

No inverse A-i(x), in a strict sense, exists when Ai(x) is non-monotonic. 
For example, the equation y = arc sin [sin x], is satisfied not only by y = x, 
but also by y = x + 2mir, and y = (2m + l)ir — x, where m is an arbitrary 
integer. For these reasons the study of the complete behavior in the 
complex plane of a transforming function B(t), for a given A(x), does not 
present the interest that might be at first expected. 

Multiple-Valued Functions. Weights. 

15. We shall now consider the case of multiple-valued real functions. 
Let y = Ai(x), y = Fi(x), y = Gi(x), • • •, be a set of single-valued func- 
tions, which taken in their totality constitute a multiple-valued function 
y = <pi(x). The expression y = <p 2 (x) may be regarded as an abridged 
way of writing the set of one-valued functions 

y = A 2 (x), y = AJLFtix)], y = A 1 [G 1 (x)] ) ■■-, 

y = F 1 [A 1 (x)] > y = F 2 (x), y = F 1 [G 1 (x)], •••, 

y = G 1 [A 1 (x)], y = G 1 [F 1 (x)] > y = G 2 (x), 

In most cases of multiple- valued functions <p, the expression, y = <p n (x), 
can have no meaning for n other than a positive integer. It sometimes 
happens that two or more of the original one-valued components of 
y — <pi(x), are commutative, e. g. perhaps Fi[Ai(x)] = Ai[Fi(x)]. In 
this case, y = Ai[Fi(x)], occurs twice in the expression y ■= (p 2 (x). This 
requires us to consider the curves as weighted, each curve having a positive 
integer weight, corresponding to the number of times that it is to be counted. 
One of the most important cases of a multiple-valued real function is 
the very special case in which <p(x) consists of two parts, one of which is 
simply y = x, and the other of which is a real monotonic, increasing, 
single-valued function, y = Ai(x). We shall obtain in this case, 

(po(x) = x =3 A (x), 

<p x (x) m {Ao(x), Ai(x)}, 

<p 2 (x) m {A (x), 2Ai(x), A 2 (x)\, 

<p 3 (x) = {A (x), SA i(x), 3A 2 (x), A 3 (x)}, 

<p n (x) = \A (x), (1)A 1 (x), (t)A 2 (x), •••, ®Ai(x), ■■■}, 

for n a positive integer. By extending the notion of weight so as to include 
negative, fractional, and irrational values of the weight, we may define 
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<p n (x) for all real values of n by the relation 

<p n (x) eb {A (x), (1)A 1 (x), ®A x (x), •••, (1)A n (x), ■■■}, 

the expressions (*) indicating as usual the coefficients in a binomial expan- 
sion. It is, furthermore, only by some such extension that the expression 
<P-i(x) can be determined or even defined. 

Matrix Interpretation. 

16. The iteration of a real function derives a great deal of added sig- 
nificance when regarded as a matrix multiplication. The matrices A n < - m) 
which we have hitherto considered have contained but a discrete set of 
elements. We shall now consider matrices of a more general type, in which 
the elements that appear in the matrix may form continuous systems. 
One such type of matrix is suggested by the Ai (m) already considered. 
The elements of the (i + l)st row of Ai (m \ are the coefficients of A^x*), 
where the coefficient of x j appears in the (j + l)st column, i and j being 
integers. If now we allow i and j to take on all real positive values, we 
shall obtain an analogous matrix, in which zeros occur everywhere except 
on the main diagonal and lines parallel to it starting with the columns at 
integral units distances from the initial element, and in the first row. 

We shall now introduce certain functions M n (x, y) of two variables, 
x and y, which functions we shall call matrices. The matrix M\(x, y) 
associated with a single valued function y = Ai(x), will be the function 
which vanishes for points not on the locus y = A\{x), but is unity on the 
locus. More generally, for a multiple-valued and arbitrarily weighted 
function y = <pi(x), we shall mean by Mi(x, y) the function whose value 
at any point (x, y) is equal to the weight of the point, where points not on 
any component of y = <pi(x) are considered as of weight zero. The matrix 
product of Mi(x, y) by itself we shall call M 2 {x, y) and in general the matrix 
product of Mn-i(x, y) by Mi(x, y) we shall denote by M n (x, y). The value 
of M 2 (x, y) for any point (x, y) is defined as the algebraic sum of the values 
of the algebraic product of Mi(x, t) by Mi(t, y), as t varies over the interval 
a, b, in accordance with the usual notion of matrix product extended in an 
obvious manner. 

If we write y = <pi{t), and t = <pi(x), and give x and y arbitrarily 
selected values, y\ and Xi respectively, the roots of yi = <pi(t), t = (h, h, t 3) 

• ■ •) are in general all distinct from the values of t = <pi(xi), t = (£/, U , t 3 ', 

• • •)• Let however the common values, if any, be denoted by </', t 2 " ', W , 

• • • . The weight of the point (</', yi) we may denote by a,-, i = 1, 2, • • • , 
and the weight of the point (xi, t"), by fr, i = 1, 2, • • •. To obtain 
y = (p 2 {x), we must eliminate t from y = <pi(t) and t = <pi(x). In particular 
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the weight of the point (xi, t/i) will be the sum of the weights given by the 
points, t\", W, t 3 ", • • • , or in other words it will be 2,-a,-/3,-. Hence the 
process of determining <p 2 is much the same as that of finding M 2 (x, y); 
in fact the matrix M 2 {x, y) is exactly the matrix M corresponding to 
y = <p 2 {x), and, in general, the matrix of <p n (x) is exactly M n (x, y). 

We may now find an entirely new interpretation to one of the formulas 
already given. Let us denote by a lx z a matrix which vanishes except for 
x = y, for which equality, it assumes arbitrary values along this diagonal. 
By Viiz) we shall denote a matrix which vanishes along the main diagonal, 
i. e. for x = y. By A 1(2) we shall denote the matrix which is the sum of 
aiiz and Vi(z). We shall regard z as a symbolic variable upon which we 
are operating with the matrices a n , Vi, A h etc. The matrix product of 
the operation a^z by itself we shall denote by au 2 z and the inverse of auz 
by (l/an)z, etc. We shall then have the formula (5) as holding also with 
this new interpretation, namely, 



Mx) = «nir.« + (1 )£/i(*) + (2)^^) 



+ 



+ 



where 



(;)>•» 



+ ■••], 



Vi(z) m Ai{z) - (;>n^_i(2) + (*)a n 2 ^i-2(2) + ■ ■ ■ ± a u VLo(z) 

and where in the present instance Vt(z) is the ith iterate of V\{z). It is 
particularly important as a formula, when V\{z) represents a curve which 
is one-valued and monotonic. The weight of Vi{x) may vary from point 
to point. When On is equal to unity along the main diagonal and the 
weight of Vi is either one or zero, Vi being the matrix of a one-valued 
monotonic function, we reobtain the formula given in § 15 for <p n (x). 

We have interpreted the iteration of a real function as finding the power 
of a matrix. It is also possible to interpret finding the power of an ordinary 
matrix as a case of iteration. For concreteness let us consider the example 
of the following normal matrix: 

n 811 

n s 12 

n 

Ai = r 2 S21 

r 2 

r 2 

n 
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Here we have for an 
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Here V 2 contains but a single element different from zero, and F 3 , V4 
vanish identically. From the formula, we obtain 

n(n — 1) 



A n m 



ri n Wi n-1 Sii 



ri"- 2 siisi 2 



r n 









nri" -1 Si2 












r 2 " 














W 2 n-1 S 2 i 

r 2 n 
















r 2 n 




which may be readily verified.* 








r 3 n 



PART THREE. 



SOME PROPERTIES OF CERTAIN RATIONAL FUNC- 
TIONS UNDER ITERATION. 

A Special Class of Rational Functions. 



17. One of the simplest and most interesting classes of functions 
y = Ai(x) to study in the complex plane, is that of the rational functions. 

* Cf . Van Vleck, "One Parameter Projective Groups . . . ," Trans. Am. Math. Soc, 
13 (1912), p. 353-386. 
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Certain rational functions with real coefficients have the property that the 
iterates of the real function are identical with the iterates considered as 
real functions. In these cases the positive and negative iterates of a real 
point are themselves real. It is a class of these functions that we shall now 
consider. 

Let us be given a rational function Ai(x) of degree r in numerator and 
denominator. For one value of x, y = Ai(x), will have but one value, 
while for y assigned arbitrarily x must have r values. If these are to be 
real whenever y is real, which we must require in order that real points 
are carried into real points by A-i, then the r zeros of Ai(x) must be real 
and distinct; since if y = yields two coincident roots, then y = c, where 
c is properly selected with respect to sign, and different from zero, will have 
two imaginary roots. Furthermore and for the same reason y = Ai(x) can 
have no horizontal tangent at a finite point, so that between any two zeros 
of y = Ai(x), there must be a pole of this function, and this continues to 
be true when we regard the range cyclically. Thus the poles and zeros 
alternate. The curve y = Ai(x) must then either increase monotonically 
between every pair of vertical asymptotes, or else decrease monotonically 

between every such pair. In either case the expression —4 — must 

have the same sign at all of the zeros of Ai(x). The straight fine y = x 
meets the curve y = A{x) in r + 1 points real or imaginary. If the curve 
decrease monotonically between asymptotes, that is, if [dA i(x)]/dx is negative 
at the zeros of Ai(x), then these r + 1 points of intersection must all be 
real and distinct, while for the curve increasing between asymptotes, r — 1 
at least must be real and distinct, these lying between successive vertical 
asymptotes, while the two remaining are either conjugate imaginary, real 
and coincident, or real and distinct. 

The sum of the residues of lj[Ai{x) — x], must be zero. We shall sup- 
pose that the points of intersection e h e 2 , • • ■, e r+h of y = x with y = Ai(x) 
are all finite. We may then obtain a relation among the slopes of the 
curve at the points e, by use of the residues of l/[Ai(x) — x\. The curve 
y = Ai(x), will be supposed to have r finite poles, and hence it must have 
a real horizontal asymptote, say, y = a. If l/[Ai(x) — x] be expanded in 
descending powers of x, its first coefficients will be given by using the 
expression of the form 



Ai(x) - x P(x) 

X Q(x) 



- 1 



~ I ~2 I ,v.3 I 

f •*/ *t/ *t/ 

'~l+5! + «l + 
xx 2 
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which determines the residue at infinity of l/[Ai(x) — x] as being equal 
to unity. If for dAi(x)/dx we write Ai'(x), the residue at e,- will be 
l/[Ai'(e,) — 1], hence we have, 

r+1 



r+1 1 

y - ~.= i 



If Ai(x) decreases between the vertical asymptotes, then Ai(x) is 
certainly negative in particular for the r + l points e,-, % = 1, 2, • • •, r + 1. 
If it is increasing between the asymptotes, then Ai'(x) is real and greater 
than unity at each of r — 1 points e, one lying between each pair of vertical 
asymptotes. If the two remaining e's are real and distinct then Ai'(x) 
will be greater than unity at one of them, and positive but less than unity 
at the other. 

In case Ai(x) is negative at the points e, there can be no question of a 
real one- valued function y = A n (x) which is a continuous function of n, 
and which leaves one of the e's invariant. For if such an A n (x) existed its 
slope at any e, i. e., A n '(x) would be equal to [Ai(x)] n at the same point, 
but this cannot be real and vary continuously with n and still be negative 
for n — 1. We shall examine now in greater detail the case when a func- 
tion y = A n (x) might exist and vary continuously with n, i. e., the case in 
which Ai'(x) is positive at the real points e. We shall confine ourselves 
also to the case in which not only r — 1, e's are real and distinct but in 
which the two remaining e's are real and all r + 1 are distinct. The cases 
in which the two remaining e's are conjugate imaginary or coincident may 
be treated in an analogous manner and present no new features of special 
interest. 

The Iteration of Certain Rational Functions. 

18. If we be given a rational function, y — Ai(x), we may, as we have 
stated before, replace x by a linear fractional transform of x, and y by the 
same linear fractional transform of y, and the transform of Ai obtained in 
this manner has for its iterates the transforms of the iterates of the original 
A i. By taking a real linear fractional operation we may arrange matters 
in such a way that r — 1, e's, e x , e 2 , • • •, e r _i, occur each between a pair of 
vertical asymptotes, while e, and e r+ i occur to the right of the right-most 
asymptote. 

We shall suppose therefore that we are given two polynomials P(x) and 
Q(x) of degree r, and we shall write Ai(x) = P(x)/Q(x). We shall suppose 
P and Q to be real, and the r zeros of P(x), z u z 2) •••,«, to be real, and the 
r poles of Ai(x), i. e., the r zeros of Q(x), pi, p 2 , • • •, p r , to be real while 
— °o < Pi < zi < ?>2 < z<i, • • • , < p r < Zr < + <*>• We shall also suppose 
that Ai(z) > for each point z. There will be r — 1 points e,, i = 1, 2, 
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• • •, r — 1, for which Ai(ei) = e,-, and such that pi < e t < p,+i, while 
Ai(e,) > 0. The two remaining roots of A\(x) = x, e r , and e r +i, we may 
denote also by /i and / 2 respectively. We shall suppose these to be real 
and distinct, and satisfying the inequality p r < /i < / 2 . Then Ai'iji) > 1, 
and < Ai'(ft) < 1. 

For n a positive integer greater than one, the function y = -A„(x) has 
many features analogous to those of y = Ai(x). It is a rational, real 
function of degree r n , with r n real zeros and r n real poles, and the poles 
and zeros alternate. The roots of A n (x) = x are r n + 1 in number of 
which r" — 1 are real and occur one each between every pair of poles, while 
the last two are real and to the right of the right-most pole, and coincide 
indeed with/i and/ 2 . The left-most pole of A n (x) is to the right of the left- 
most pole of A n -\{x) but lies also in the interval p h p 2 . The right-most pole 
of A n (x) lies to the right of the right-most pole of A n _i(x) but also in the 
interval p r fi. The function y = A n (x) increases monotonically between 
every pair of poles, while A n '(fi) > 1 and < A n '(f 2 ) < 1. If ai denotes 
the asymptotic value of y = A x (x), then the asymptotic value, a„, of A n (x), 
is such that a n < «i. 

For n positive but nonintegral, the functions A n {x) which we can 
determine in the neighborhood of /i and of / 2 coincide, in the sense that the 
series valid at /i when analytically extended must be valid at fo, and must 
coincide with the series obtained at / 2 , for no singularity can occur between 
these two points. It might be at first supposed that by further continua- 
tion we should obtain a function with essential singularities, perhaps, on 
the real axis, but none the less one-valued in the complex plane. This 
cannot however be the case if n be allowed to vary continuously, as we 
shall proceed to show by the use of the Riemann surface for Ai(x). That 
y — A n (x) is not one-valued when n = — 1, and r (the degree) is greater 
than one, is obvious. 

Let us mark in the complex plane of x, the branch points of y = A- X (x), 
and join these by cuts, to obtain a region R , which is simply connected and 
includes all points of the complex plane not themselves on the cuts. Since 
/i is not itself a branch point we shall suppose that no cut actually passes 
through /i, and likewise for / 2 . Now the interior of a small circle about 
U is carried by the transformation y = A-i(x) into the interior of a curve 
approximately circular and concentric with the given circle, since f% is 
carried into itself, and A_/(/ 2 ) is real and greater then unity. So that any 
small region, X, inclosing / 2 is carried into a region Y entirely including X 
within its interior. Now A-i(x) is denned as one-valued, either finite or 
infinite, for every x in R , so that as X is successively extended Y suffers 
a like extension and always so as to include X, as long as X is within Ro- 
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As X comes to coincide with R , Y is, of course, no longer a non-overlapping 
region but coincides with the Riemann surface for y = A-i(x), except for 
certain cuts, a region which we shall call R-i. We shall define itL 2 as the 
multiple-sheeted region that Y becomes when X is taken as R-i, and in 
general for X = R- H we shall have Y = R-( n+ i). On the other hand, 
for Y = R 0) X becomes a region included within R and only partially 
covering the complex plane. This we shall call R\. We shall define anal- 
ogously jR 2 , #3, • • • • Now we have a series of regions, 

• • •, .R3, Ri, R\, Ro, R-i, R—2, R-3, • • • 

each included in the following, and building up successively an infinitely 
sheeted surface, which we shall call R. 

By the discussion made in the first part of this paper, we may obtain 
the coefficients of a power series E(x) which represents an analytic function, 
such that E[A-i(x)] = (l/a n )E(x). If we write t = E(x), then the trans- 
formation y = A-i(x) corresponds in the i-plane to the transformation 
t' = (l/a n )t where t' = E[A-i(x)]. Thus the mapping of R n upon 2?„_i 
has for its transform in the £-plane, a mere magnification leaving the origin 
invariant. The straight lines through the origin in the i-plane are carried 
into themselves by the magnification, and are the path curves of the trans- 
formation in the i-plane for t' = (l/au n )t. The corresponding path curves 
in the x-plane, are the path curves for A- n {x), as n varies over real values, 
and are curves starting at / 2 and passing, in their totality, through every 
point of R. For t = E(x), as t varies over the i-plane, x varies over R, 
so that R is the Riemann surface for t = E(x) . Now let us take a number x. 
There will be one point in each sheet of the infinitely-sheeted surface R, 
having x for its coordinate. Through each of these points there will be a 
path curve terminating at/ 2 . To determine A n (x), for a given real, positive 
n, we must follow along a path curve from x toward / 2 a distance correspond- 
ing to n in the scale that is determined along the path curve by its very 
nature. The infinite number of points which we shall obtain by starting 
with the same x but in different sheets, can only be superposed for all 
choices of x and n when the path curves in the different sheets of R are all 
superposed. This is not, however, the case unless the rational function be 
linear fractional, since otherwise y = A-i(x) would carry a single x into 
but a single y. Hence A n (x) is not a one- valued function of x as n varies 
over real positive values. 

The function x = E-i(t) is a one- valued meromorphic function of t 
with infinity as an essential singularity, while E{x) itself is an infinitely 
many-valued function with R for its Riemann surface. 

The same sort of discussion will show that if A (x) be an analytic function 
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defined with its inverse, all over the plane and Ai(x) is one- valued, then 
while A n (x) cannot be one-valued for general values of n, still the inverse 
E-i of the transforming function E is in general one-valued and mero- 
morphic over the complex plane, provided that there is one point of / used 
in defining E for which A X (J) = /, and | A^'if) | = A/(/) < 1. "We shall 
not investigate less stringent conditions at this time. 

Pbinceton, N. J., 
February, 1915. 



